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Abstract. First, we briefly outline some aspects of the starting project to design non-empirical 
energy functionals based on low-momentum vacuum interactions and many-body perturbation 
theory. Second, we present results obtained within an approximation of such a scheme where the 
pairing part of the energy density functional is constructed at first order in the nuclear plus Coulomb 
two-body interaction. We discuss in detail the physics of the odd-even mass staggering and the 
necessity to compute actual odd-even mass differences to analyze it meaningfully. 
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INTRODUCTION AND ELEMENTS OF FORMALISM 

m 

Like-particle pairing is an essential ingredient of nuclear-structure models, in particular 
regarding the description of exotic nuclei yj]. Also, superfluidity plays a key role in neu- 
tron stars, e.g. it impacts post-glitch timing observations [0] or their cooling history (H]. 
Within a single-reference (SR) implementation of the energy density functional (EDF) 
Q\ ', formalism [j4|], pairing is incorporated through the breaking of the U(l) symmetry 

associated with particle-number conservation. As a result, the binding energy S'sr of 
the many-body system is postulated to be a functional of both the one-body density 
matrix pji = (<J>|cJc/|<E>) and the pairing tensor Kji = (4>|c,c/|4>), the dependence on the 
latter being allowed by the use of an auxiliary product state of reference |4>) that mixes 
particle numbers (of given parity) |0]. 

Modern empirical parameterizations of existing EDFs, e.g. Skyrme or Gogny, provide 
a fair description of bulk and certain spectroscopic properties of known nuclei Q4Q _ On the 
other hand, they lack predictive power away from known data and a true spectroscopic 
quality, in particular regarding the part that drives superfluidity. As a result, several 
groups currently work on empirically improving the analytical form and the fitting of 
functionals, e.g. see Refs. [@,0] for recent attempts to pin down the isovector content of 
purely local pairing functionals. 

Along with improving the phenomenology at play, the quest for predictive EDFs starts 
to benefit from a complementary approach [8] that does not primarily rely on fitting 
known data but that roots the analytical form of the functional and the value of its cou- 
plings into underlying low-momentum two- and three-nucleon (NN and NNN) interac- 



tions [SOU through the application of many-body perturbation theor>Q (MBPT) iflill . 
The overall goal of such a project is (i) to bridge with ab-initio many-body techniques 
applicable to light nuclei, (ii) calculate properties of heavy/complex nuclei from basic 
vacuum interactions and (iii) perform controlled calculations with theoretical error bars. 
First results following such a route are currently being reported lfl5ll . It is an objective 
of the present contribution to expose results of such an effort to build the pairing part of 
the EDF non-empirically H O, EE EE . 

We propose to write the energy functional at a given order in (Goldstone) MBPT 
under a generic form that is convenient to bridge with existing phenomenological EDFs 

<?SR[{Pij},{Kij},{Kijh{Ek}} = J>P* (1) 
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where all dependencies on p and k*k have been made explicit. The effective vertices 
v P j kl , vfjki- ■ ■ thus introduced are expressed in terms of the vacuum two-, three-,. . . body 
interactions and on quasi-particle energies that are to be determined self-consistently 
through a chosen procedure. More precisely, a term of given power in p and/or k*k in 
Eq.[T]receives contributions from different perturbative orders and/or many-body forces. 
To exemplify this, we can write the vertices arising at second order in the NN interaction 
v^,in a perturbation theory that does not account for pairing explicitly^ 
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where denotes single-particle energies to be determined self-consistently. The EDF 
form of Eq. [T]may naively suggests that it results from the average value, in the unper- 
turbed vacuum, of an (hypothetical) effective Hamilton operator containing two-body 
(second line), three-body (third line),. . .pieces. However, Eq. [2l that provides micro- 
scopic expressions for the matrix elements of v^ z , vf^, ^ k P mn - ■ ■ , demonstrates that 
re-extracting an (effective) Hamilton operator from the energy density has no founda- 
tion^ and can at best be the result of approximations. 



Infinite resummation of certain categories of diagrams and/or a redefinition of the unperturbed vacuum 
|<I>} are always possible. Switching from conventional hard-core potentials to low-momentum interactions 
is essential to make a perturbative approach viable, e.g. second-order calculations performed in terms of 
low-momentum interactions provide satisfactory results for bulk correlations ifTTl [l2l [l3ll . 

2 The single-particle basis solution of Eq.[3]also diagonalizes the density matrix p of |<J>) in this case. 

3 Note for instance that symmetry properties of v^L, v^ mn and vf^^ w under the exchange of fermionic 
indices are not as expected from two-, three- and four-body operators. 



Forms as given by Eq. Q] are known as orbital-dependent energy functionals 112 lh in 
electronic systems density functional theory (DFT), with the important subtlety that DFT 
implies that quasi-particle wave functions (C4,\4) and quasi-particle energies are 
generated through the variationally optimum local one-body potential, i.e. the optimal 
effective potential (OEP) I122ll . We do not insist on that here to rely on a framework 
that embraces empirical Gogny functionals whose associated one-body fields are non- 
local0. Of course, none of the existing empirical functionals do depend on quasi-particle 
energies E^. It remains to be seen in the future whether such an extension is necessary 
and tractable. 

An alternative to OEP that is closer to what is currently done with empirical EDFs 
consists of determining quasi-particle wave-functions and energies through the mini- 
mization of Ssr with respect to independent matrix elements of p and K, under the 
constraint to have given neutron and proton numbers in average. This leads to solving 
Hartree-Fock-Bogoliubov-like (HFB) [5] equations 

( ^a* -ff+k )\vl) =Ek (vkJ- (3) 

The one-body field h that drives the correlated single-particle motion and the shell 
structure, as well as the field A that drives superfluidity, are defined as 

°Pn u ° K ij 1 u 

through which two effective vertices v ph and v pp are introduced that can be expressed in 
terms of v^, vf^. . . , i.e. they themselves possess a diagrammatic expansion in terms 

ofv^andv™. 

Our immediate focus is on the pairing part of the EDF Beyond enhancing its predic- 
tive power, our aim is to understand better the microscopic processes that build superflu- 
idity in finite nuclei. Typical questions relate to (i) the contribution from the direct NN 
and NNN interactions, its breaking down in partial waves (essentially ^o, 3 ^i, ^2 in 
decreasing order of expected importance), as well as (ii) the role of higher-order effects 
associated with the coupling to (collective) fluctuations. 

To answer the first of these two questions, our current target is to perform reliable 
finite-nuclei calculations at first order in low-momentum NN and NNN interactions gen- 
erated through renormalization group techniques [^,113]. The upper row of Tab. Q] shows 
the corresponding diagrammatic for the one-body fields, omitting for simplicity contri- 
butions from the NNN interaction. As an intermediate step, we present here approximate 
results such that v ph and the part of the EDF that depends only on p is empirically pro- 
vided by the SLy4 Skyrme parametrization 112411 . while pairing vertices v KK = v pp are 



Using a non-local pairing field, as in the present work, renormalizes from the outset the ultraviolet 
divergence that arises when using a (quasi-)local pairing field. 

5 Not only the present scheme does not insist on obtaining a local potential but also partial derivatives 
with respect to quasi-particle energies are omitted. The so-called Rrieger-Li-Iafrate approximation to the 
OPE 1 23] also omits such functional derivatives. 



TABLE 1. Perturbative expansion scheme to first (up) and second 
(down) order. The dashed line denotes the free-space NN interaction. 
Diagrams with more than one anomalous propagator are not shown. 
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computed at first order in the Coulomb plus nucleatO NN interaction [16, 17, 18, 191. 
Such an EDF contains neither energy dependencies nor p-K cross terms. Only the dom- 
inant x Sq partial-wave of the NN is included in v KK whereas the effect of 3 Pi and 1 Z>2 is 
discussed in Ref. Ir26ll . Also, the first-order contribution of the NNN interaction to v pKK 
and v pp will be reported on in Ref. [20]. Note that for such a calculation to be a de- 
cent approximation of the targeted first-order one, at least as for extracting pairing gaps, 
it is crucial that the empirical Skyrme parametrization that drives the underlying shell 
structure is characterized by an isoscalar effective k-mass ps 0.7 m at saturation den- 
sity [|19|l . Of course, we eventually aim at calculating h at lowest-order in both the NN 
and the NNN low-momentum interactions, possibly making use of the density matrix 
expansion [27[|28]. Eventually, higher-order contributions are left out for future workl! 



EXPERIMENTAL VERSUS THEORETICAL PAIRING GAPS 

We limit ourselves to discussing the odd-even mass staggering (OEMS) whereas other 
observables are reported on in Ref. [20]. The OEMS is dominated by the deficit of 
binding energy of the unpaired nucleon in odd nuclei, i.e. the "pairing gap". In the SR- 
EDF formalism, such a staggering relates to the description of odd nuclei through the 
excitation of a quasi-particle on top of an even-number parity vacuum. Experimentally, 

the OEMS is extracted through n -points difference-mass formulae A^(iV/Z) S 



However, because of the technical difficulty to compute odd nuclei, data are often 
compared to purely theoretical pairing gaps extracted from the calculation of just one 
even-even nucleus. One such theoretical pairing gap is the Lowest Canonical State gap 
Alcs> defined as the diagonal matrix element of the pairing field A in the canonical 
single-particle state whose associated quasi-particle energy is the lowest iflvL fl~8h . 
The difficulty with such comparisons is not only that (i) finite-difference mass formu 



lae are contaminated by contributions other than the targeted "pairing gap" 00113 11 132 , 



331133] but also that (ii) the "pairing gap" that makes the actual OEMS is itself an aver- 



age of Alcs extracted from the even-even and (blocked) odd-even nuclei involved in the 



finite-difference mass formula |33L I34f |. As a result, comparisons based on theoretical 



gaps extracted from one even-even nucleus can only be of semi-quantitative character, 



6 We use the low-momentum NN interaction Vi ow k Hi built from the Argonne v\g NN potential ll25ll at a 
renormalization cut-off A = 2.5 fm _1 . 

7 See the lower row of Tab. Q]for the second-order contributions to the one-body fields. 
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FIGURE 1. Experimental three-point mass differences (crosses) for neutrons along the tin isotopic 
chain versus several theoretical measures of the OEMS: A£ cs (even) (dashed-line), A£ cs (odd) (dashed- 

(3) 

dotted line) and A n (N) (full line). Theoretical three-point mass differences are also shown for odd- 
even nuclei computed using an even-number-parity vacuum as a reference state, i.e. a HFB state without 
any quasi-particle blocking as if odd-even nuclei had the same structure as even-even ones (dotted 
line) IM HI. 



which is often fine as empirical pairing functionals are not yet targeting a nucleus by 
nucleus agreement with experiment. However, aiming at such a level of agreement in 
the (distant?) future and at doing so in a non-empirical fashion requires the comparison 
of apples with apples, i.e. to compare theoretical and experimental odd-even mass differ- 
ences. At the price of requiring a good understanding of the different contributions to the 
OEMS [|33L l34f| . doing so allows more fruitful comparisons between theory and exper- 
iment, e.g. to analyze the interplay between pairing and the underlying shell structure. 
This is what we wish to briefly exemplify in the present contribution. 



RESULTS 



In Refs. iQjLGji], experimental (odd) were compared to A? cs (even). Main results 



were that neutron and proton pairing gaps computed from the Skyrme plus non-empirical 
pairing energy functional were close to data for a large set of semi-magic light-, medium- 
and heavy-mass nuclei. Implications of such results were also discussed. Here, we wish 
to analyze the qualitative modifications brought about by comparing directly theoretical 
and experimental three-point mass differences. To do so, we computed odd-even nuclei 
through the self-consistent blocking procedure performed within the filling approxima- 
tion [|35L[36|1 . Results for neutron gaps along the tin isotopic chain are reported in FigCQ 
To analyze meaningfully the OEMS [13411 . the ground state of odd nuclei is best 
understood as a quasi-particle excitation on top of an even-number parity vacuum that 



shares the structure of even-even nuclei but that has the odd number of particles on 
average. In this way, the quasi-particle excitation is performed at (almost) constant 
particle number [33]. The even-number parity vacuum provides the smooth part of 
the energy while the quasi-particle excitation, that is dominated by the static "pairing 
gap", generates the relative lack of binding of odd nuclei without which no interesting 
physics would be extracted from odd-even mass differences. The dotted line in Fig. Q] 

(3) 

shows the contribution of the smooth part of the energy to A„ (N), i.e. when odd- 
even isotopes are described as if they had the structure of even-even ones. One sees 
that such a contribution, which reflects the curvature of the smooth part of the energy, 
oscillates symmetrically around zero and accounts exactly for the odd-even oscillation 

(3) 

of A„ (N). This demonstrates that contributions other than the targeted "pairing gap" 

(3) J — li 

contaminate A„ '(N) in an opposite way for odd and even N 113411 - which contradicts the 

usual belief r[30h that A^ (odd) is free from such contaminations^. 

In Fig. [H the comparison between experimental data for odd N and AL CS (even) 
(dashed line) recalls the results of Refs. fTTvL [l8ll and sets the stage for what comes 
next. Those two curves are consistent with each other, with a slight overestimation 
(underestimation) of the data at mid-shell (just below and above the N = 82 shell 
closure). Such a situation is representative of the results obtained along other semi-magic 
isotopic and isotonic chains. Still, certain features that are visible in the data, i.e. (i) the 
lowering around N = 65, (ii) the flat trend as one approaches the N = 82 shell closure 
and (iii) the finite jump from A = 81 to A = 83, are not reproduced by the bell- shaped 
curve provided by A^ cs (even). At best, one can talk of an overall semi-quantitative 
agreement and wonder whether the remaining discrepancies are due to limitations of (a) 
the pairing part of the EDF, (b) the use of A^ cs (even) as a measure of the OEMS and/or 
(c) SR calculations that miss dynamical correlations associated with particle number 
restoration and collective pairing vibrations that are of importance in the weak pairing 
regime, e.g. near shell closures. 

The full-fledged comparison of experimental (stars) and theoretical (full line) three- 
point mass differences is also provided in Fig. [TJ The most striking feature is the 
ability of the calculation to grasp quantitatively the three non-trivial features seen in 
the data and outlined in the previous paragraph. As a result, one goes from a semi- 
quantitative agreement with experiment across the major shell using A^ cs (even) to the 
ability to compare on a nucleus by nucleus basis. In particular, there were hints that the 
lowering of the pairing gaps around N = 65 could be partly due to dynamical pairing 
fluctuations Il37m . Here, such a feature is well reproduced at the SR-HFB level. It will be 
of interest to study whether using the non-empirical pairing functional computed from 
the finite-range and non-local Vi ow k interaction is essential to obtain such a pattern or if 
it is entirely driven by the interplay with the underlying shell structure, independently of 
the detailed characteristics of the pairing functional employed. 



It was suggested in Ref. 113411 to use A„ (odd) as a measure of the sole "pairing gap" because the 
contribution from the time-odd reversal symmetry breaking, not discussed in the present paper, possibly 
cancels out the contribution from the smooth part of the energy in this case. This is however subject to 
revision due to the current lack of knowledge regarding time-odd terms in the nuclear EDF. 



A similar situation occurs regarding the behavior of the OEMS towards and across 
the N = S2 shell closure. One sees from A^ cs that static pairing correlations collapse in 

the immediate vicinity of N = 82, i.e. in 13 l,l32,l33g n ^ w Yiile fa e experimental A„(N) 
sustains a non-zero value down to N = 81 and 830. Surprisingly enough, when going 
from A^ cs to theoretical three-point mass differences, the experimental trend is well 
captured down to N = 81 and across the N = 82 where the OEMS jumps by 350 keV, 

An (83) being the last piece of available data. Due to the collapse of A^ cs close to 

(3) 

the shell closure, it is usually stated that A„ is dominated by other contributions than 
static pairing in this regime, i.e. by dynamical pairing fluctuations and contributions 
associated with the discreteness of the underlying shell structure. Regarding the former, 
we just saw that a SR calculation omitting entirely dynamical pairing fluctuations^ can 
account for the data. Regarding the latter, it is to be noted that (i) in the (hypothetical) 

(3) 

zero-pairing limit, and still assuming spherical symmetry, A„ (N) is zero from N = 70 
to N = 81 as one fills the highly degenerate h n / 2 snei l and that (ii) the regularly 
oscillating contribution of the smooth part of the energy seen in Fig. Q] demonstrates 
that the structure of odd nuclei is still best understood, down to 131 Sn and 133 Sn, as 
a quasi-particle excitation on top of a statically paired even-number parity vacuum. 

(3) 

Eventually, the energy of the quasi-particle excitation that builds A„ (81, 83) and leads 
to the unpaired blocked state is dominated by pairing correlations, i.e. it would be zero in 
the zero-pairing limit. Although dynamical pairing correlations are likely to renormalize 
the OEMS, the present results implies that odd-even mass differences might be less 
impacted by such correlations than other observables in the vicinity of shell closures. 

CONCLUSIONS 

We discuss pairing gaps obtained in tin isotopes using an energy density functional 
whose pairing part is constructed at first order in the nuclear plus Coulomb interac- 
tion. Only the (dominant) [ Sq partial wave of the two-nucleon force is incorporated 
whereas the contributions from 3 P\ and l D 2 ||2~6ll . as well as from the three-nucleon 
interaction |20j], will be reported on soon. Most importantly, we discuss in detail the 
physics of the odd-even mass staggering and the necessity to compute actual odd-even 
mass differences to analyze it meaningfully and compare with data on a nucleus-by- 
nucleus basis. In particular, an excellent description of the odd-even mass staggering is 
obtained in the vicinity of magic shell closures prior to incorporating dynamical pairing 
correlations associated with particle number restoration and pairing vibrations. 



9 One must remove A n (82) from the analysis as it measures the N = 82 shell gap rather than static 

(3) 

pairing correlations. Contrarily, A„ (81,83) are not influenced by the N — %2 shell gap. 

10 The Lipkin-Nogami procedure is not used in the present calculation. 
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